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Propagation of gravitational and acoustic plane waves in a flat universe filled with a general 



relativistic, homogeneous and isotropic, spatially flat continuum is studied. The continuum is 
, described by analogues of nonrelativistic characteristics, namely energy per particle, pressure 

and Lame coefficients, and considered in the comoving proper-time gauge. For all modes with 



the given wave covector, differential equations governing the time dependence of the amplitudes 



are derived. In particular, longitudinal acoustic waves are described, in analogy with the non- 
relativistic theory, by two coupled first-order equations. As an example, plane waves in a stiff 

— «. . ultrarigid continuum are considered. 

O . 

cr 

1 Introduction 

Elasticity, as a theory of matter in form of an (ideal, dissipationless) elastic continuum, 
has many applications in everyday life; but in areas where relativistic effects are important, 
applications are much more sparse. Nevertheless, the presence of a solid crust in neutron stars 
has motivated the development of relasticity (short name for relativistic elasticity, from p]). 
Among early works on the topic, [2] appears to be the most successful. 

Since then, fundamental concepts of relasticity are known, but the terminology, notation 
and formalism are far from being stabilized. Different authors start from different notions of 
nonrelativistic elasticity, and use different coordinates, depending on what they find most natural 
and/or appropriate for the problem they study. This paper is no exception. We deal with a 
particular problem, thus we consider particular notions to be basic, develop particular formalism 
and perform computations in particular coordinates. 

The motivation of this paper is that, as long as elastic continuum can be included into 
the framework of general relativity, we should be able to describe how nonrelativistic elastic 
phenomena and Einstein's gravitational phenomena affect each other. 
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As a nonrelativistic effect to start with, we have chosen the propagation of acoustic waves; 
specifically, its simplest case, the propagation of plane waves in a homogeneous and isotropic 
continuum. For details, see [3]. 

In general relativity, the effects described by the theory also include plane waves in a homo- 
geneous and isotropic background, namely gravitational waves in a flat expanding universe. 

The first goal of this paper is to derive equations of propagation of both acoustic and grav- 
itational (weak) waves in a homogeneous and isotropic, spatially flat, but otherwise completely 
general continuum. Since less specific equations for general weak waves were formulated, and 
the correctness of their nonrelativistic limit was proven in 0], we do not address the problem of 
nonrelativistic limit here. 

The second goal of this paper is to show how to use the theory of relasticity when solving 
a particular problem, if simplicity and clarity are considered more important than elegance and 
universality. This goal includes introducing such formalism, and using such coordinates, that the 
result will be obtained in the easiest possible way, even if that formalism and those coordinates 
will not prove useful in other relasticity problems. 

This paper is motivated neither by astrophysics nor by cosmology, although its results can be 
applied within the latter, when considering a solid dark matter component as in [5]. Therefore, 
when describing our results, we remind the reader also of the terminology used in the theory 
of cosmological perturbations, see for example [6j. Our notations in relasticity are chosen to be 
compatible with [7], and our notations in differential geometry to be compatible with [8]. 

The paper is organized as follows. In section [2] we introduce basic concepts of relasticity. In 
section [3] we develop a theory of plane waves as linear perturbations imposed on the background 
solution. We derive equations for the background as well as for the perturbations, identify modes 
that correspond to gravitational and acoustic (transversal and longitudinal) waves, and propose 
a gauge invariant description of longitudinal acoustic waves. For better understanding of our 
formalism, in section H] we write down equations for plane waves in a particular continuum, 
namely the stiff ultrarigid continuum introduced in [9j. In section [5] we discuss the results. 

2 Relasticity 

Relasticity is a theory that describes elastic continua within Einstein's theory of gravity. In 
the simplest case there is just one continuum, and it extends over all spacetime. The spacetime 
is, of course, a 4-dimensional manifold, equipped with a metric g that is regarded as a dynamical 
field satisfying the Einstein equations. 

Properties of the continuum are defined on the material manifold, that is, on a 3-dimensional 
manifold whose points represent particles of which the continuum consists. The way the con- 
tinuum is distributed in spacetime is described by the mapping / that maps spacetime on the 
material manifold, telling us which particle is in the particular worldpoint. It is called material 
mapping and considered a dynamical field on the spacetime, too. Contravariant metric in any 
worldpoint can be push-forwarded by the material mapping, providing us with a symmetric, 
double contravariant tensor on the material manifold. This tensor can be interpreted as the 
strain tensor. 

The only characteristic of continuum that influences its dynamics is the state function g. It 
is a scalar function with two arguments, the point of the material manifold and the strain tensor 
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at that point. By definition, it is the energy density in the rest frame of the given particle, in 
the continuum whose deformation is described by the given strain tensor. 

The dynamics is determined by the action S = J {—p+R)uu g , where R is the scalar curvature, 
p is the energy density given by the state function and uj g is the 4-form of volume defined by 
the metric g (for details, see [8]). We use a system of units in which the factor in front of R 
involving gravitational constant is absent. In these units, the Einstein equations assume the form 

The principle of extremal action provides us with the equations of motion. Variation of the 
action with respect to the metric leads to the Einstein equations with an explicit formula for the 
energy-momentum tensor (to be given in the next section) . Variation of the action with respect 
to the material mapping yields equations of motion of the continuum (also written in the next 
section). As a matter of fact, Einstein equations together with the Bianchi identity usually carry 
enough information to describe the dynamics completely; nevertheless, it is useful to write down 
the equations governing the dynamics of the source explicitly. 

The action, and therefore also equations of motion, is covariant with respect to the change 
of coordinates in the spacetime, and also in the material manifold. This is quite a large gauge 
freedom. This freedom, however, can be used for fixing coordinate-dependent expressions in such 
a way that they become simpler (although less general). In this paper, we fix the gauge almost 
completely in a way we find most suitable for our problem. 



3 Plane waves 

We will investigate the propagation of weak waves, which can be regarded as linear pertur- 
bations to a certain background solution. 



3.1 Description of the continuum; the background solution 

For simplicity, we choose the elastic continuum to be homogeneous, isotropic and flat, which 
means that there exists a set of coordinates X A in the material manifold (where the material 
index A, as well as any other upper case Latin index, runs from 1 to 3) such that the state 
function is invariant with respect to Euclidean translations and rotations of X A . 

The strain tensor has coordinate components (with respect to X A ) which can be treated 
as coordinates on the space of double contravariant tensors on the material manifold. These 
components will be denoted by H AB . The true strain tensors are always symmetric, but for the 
purposes of differentiation we can extend the state function also to nonsymmetric tensors. This 
gives correct results if we let g depend only on the symmetric part of H AB , or if we symmetrize, 
when differentiating with respect to H AB more than once, the derivatives with respect to all 
H AB, s but the last (and of course, insert symmetric values into the resulting expressions). We 
will adopt the latter procedure. In general, g can be differentiated also with respect to X A . 
However, the invariance with respect to translations implies dg/dX A = 0, therefore only the 
derivatives of g with respect to (the nine components of) H AB enter the equations. 

Having fixed the coordinates X A on the material manifold, we can express the material 
mapping / through the coordinate component functions f A . Together with arbitrary coordinates 
x 11 in the spacetime (where the spacetime index fi, as well as any other Greek index, runs from 
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to 3; we will also use lower case Latin letters to denote space indices running from 1 to 3), 
we now have everything ready for the variation of the action with respect to f A . This variation 
leads to the equations of motion of a homogeneous isotropic flat continuum 

( _ d 2 g pA VK c D \ B _ 

\dH AB dH AC dH BL>9 9 I xI k ) ^ ~ ' 

By varying the material part of the action with respect to g^ u , we obtain a formula for the 
energy-momentum tensor, 

-^V = ^ gjjAB f,i*f,v ~ 89 nw (2) 

Having completed a general description of the continuum under consideration, we can pro- 
ceed to further simplifications. The gauge freedom in spacetime coordinates can be reduced by 
requiring that f A (x) = x A {x A are comoving coordinates), after which there remains only one 
free coordinate a; . Further restriction will be goo = — 1, which implies that x° equals the proper 
time t of the particles of the continuum. The only freedom left consists in the choice of the 
hypersurface a; = 0. 

Now we can specify the background. We assume that it is, just like the continuum, ho- 
mogeneous, isotropic and flat, that is, equipped with a flat Robertson- Walker metric g — 
—dt <g> dt + a 2 (dx 1 <g> dx 1 + dx 2 (g> dx 2 + dx 3 (g> (fa 3 ), where the scale parameter a depends only 
on t. The symmetry allowed us to choose the zero time hypersurface to be orthogonal to the 
worldlines of all particles. Instead of the scale parameter, it turns out to be more convenient to 
use its logarithm z = In a. Then we have goo = — 1, .9oi = 0, g%j = e 2z 5ij, and we can compute 
the Einstein tensor to obtain 2G*oo = 6i 2 , 2Goi = 0, 2Gy = e 2z (— 4z — 6z 2 )Sij, where the overdot 
means differentiation with respect to t. The energy-momentum tensor simplifies to Too = Q, 
Toi = 0, Tij = 2dg/dH l: > — e 2z gSij. Equations of motion of the continuum hold identically. 

The strain tensor of the background solution is H AB = e~ 2z 5 AB , and the invariance of g 
with respect to rotations ensures that the derivatives of g with respect to H AB , evaluated for the 
background solution, are proportional to Sab (the first derivative) or equal to linear combinations 
of the products of Sab (the higher derivatives). Here, S is in fact the Euclidean metric of the 
material manifold, but in our coordinates we can use notation reminding us that Sab and 5 AB 
are unit matrices. This enables us to employ an extended version of summation convention, 
for example g AA — g AB SAB and gAA — gAsS AB . Since we have, in fact, identified the spatial 
coordinates with the material coordinates, we can pull-back 5 to obtain the unit matrices Sij 
and <5 y . The first matrix was already used in the formulas for 2Gy and . 

As the next step, let us introduce material characteristics to appear in the equations for 
the background, as well as in the equations for perturbations. First, Einstein's cosmological 
constant can be treated as the vacuum energy density po contributing a constant term to g. 
After separating out this term, we can write the remaining energy density as a product of the 
particle density n and the energy per particle e. In this way we obtain the general formula 
q = po + ne, which can be written for the background metric as g = po + e~ 3z £- Note that the 
separation of the vacuum energy density is in a sense artificial. We could choose any other value 
of po and change e in such a way that g will be the same. Nevertheless, the introduction of po is 
useful since it provides us with a simpler formula for the quantity e than we would obtain if we 
put simply g — ne. 
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In the nonrelativistic theory, when a homogeneous and isotropic continuum is deformed in a 
homogeneous and isotropic way, its stress tensor is given by one scalar quantity - the pressure 
a; and when it is relaxed, its elastic properties are determined by two scalar quantities - Lame 
coefficients fi and A (see, for example, [3]). Analogically, for a general relativistic homogeneous 
and isotropic continuum we define the scalar a by the relation 

2 g^AB = e 2Z °tAB, 



and the scalars \x and A by the relation 

8 2 e d 2 e 



dH AB 8H CD dH AB dH DC 



e z [XSabScd + p{$ac<>bd + SadSbc)] ■ 



With all this new scalars, the derivatives of g evaluated for the background metric can be written 

as 

2 ^AB=^ Z ^ + ^AB (3) 

and 

dHABQHCD + QH^h^ J = eZ ^ + ^ + ^AB6cD + ^-s)(S A cS BD + S A n8 B c)]. (4) 

The scalars e, a, fj, and A as functions of z are not independent, because the derivative with 
respect to z corresponds to the derivative with respect to H AB . The relations between these 
scalars are de/dz = —3a and da/dz = — (2a + 3A + 2/i). 

Now we are ready to write down the equations for the background. Too becomes po + e~ 3z e, 
which gives 6i 2 = po + e~ 3z e, and Ty becomes e~ z aSij — e 2z poSij, which gives, when combined 
with the equation for i 2 , —A'z = e~ 3z (e + a). If we take into account the identity for de/dz, we 
can see that the latter equation is nothing but the time derivative of the former equation. (At 
least, this holds if £ ^ 0, which we will assume in what follows in order to avoid singularities in 
our expressions. By the same reason, we will assume e + it^O.) Thus, the background is a flat 
Robertson- Walker metric satisfying 6i 2 — po + e~ 3z e and de/dz = —3a; in other words, a flat 
universe evolving as if it contained, in addition to the dark energy with the density qq, an ideal 
fluid with the density e~ 3z e and pressure e~ 3z a. 



3.2 Linear perturbations 

Having fixed the background solution we can start to inspect linear perturbations. Denote 
perturbations of the fields / and g by Sf and Sg. Of course, not every perturbation is permitted 
by the dynamics, but on the other hand, some perturbations correspond merely to coordinate 
changes and are not of physical interest. 

First, we reduce the gauge degrees of freedom in a way analogical to how we fixed the coor- 
dinates in the background solution. We require that for the perturbed solutions the spacetime 
coordinates are comoving, too, which implies Sf A — 0. We want also the x° coordinate to 
measure the proper time of the particles (x — t), which implies 6 goo — 0. This leaves us with 
the freedom of choosing the zero time hypersurface. But, unlike the background, the perturbed 
metric is generally not symmetric enough to allow us to pick a certain nice hypersurface (like the 
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hypersurface orthogonal to the worldlines of all particles). Therefore we proceed further with 
this last ambiguity left. 

Before writing down the linearized dynamical equations, we can take advantage of the fact 
that the linearized equations are (of course) linear. Therefore their solutions form a linear space 
(they satisfy the superposition principle), and can be complexified. The complex solutions have 
no direct physical meaning, but they can simplify the upcoming expressions a bit. 

Since the solutions form a linear space, it is sufficient to investigate some basis of them to 
obtain the complete picture. A nice ansatz for such basis follows from the translation symmetry 
of the continuum. This symmetry ensures that there are solutions of the form of (complex) 
plane waves, that is, depending on the space coordinates only through the multiplication factor 
cxp(ikAX A ), where k is a constant real wave covector on the material manifold. When dealing 
with just one plane wave, we can always utilize the rotational invariance of the background 
solution to rotate coordinates in such a way that x 1 becomes parallel to k, which means that 
the coordinate components of k become (k, 0, 0) (the first A; is a covector, the second k is a real 
number) . 

If we find all plane wave solutions, we can use them to construct any solution of the form of 
a (space-) tempered distribution. This means that we actually do not have a basis, because we 
miss the solutions that grow faster than polynomially with the space distance, as is the case for 
surface waves. We duly restrict our attention to the space-tempered distributions, referring to 
the fact that because of causality such solutions are sufficient for describing linearized dynamics 
in a bounded region of spacetime. 

Plane waves provide a significant simplification of our problem since the space derivatives 
acting on the perturbations become just powers of k. In further computations, some typical 
combinations of derivatives with respect to time or space coordinates occur repeatedly, so it is 
convenient not to use Sg directly, but to replace it by more complicated expressions in order to 
obtain simpler looking results. For a plane wave with the wave covector k already rotated into 
the x 1 direction we define a tensorial-looking, time-dependent quantity h by putting 

Sgoi(x) = -ikh Qt (t)e lkxl , 

5g ij (x) = e 2z Wh ij (t)e ikx \ 

Now, everything is ready for a straightforward but lengthy variation of dynamical equations. 
We will write down just the final expressions. Variation of the energy-momentum tensor gives 

5T 00 = -^e~ 3z (s + a)h k k, 

ST 0i = ikp h 0l + ike~ 3z sh Ql , 
STtj = -ie~ z (A + <j)5 i: jh kk - e^ z ph lj - e 2z pohij; 
variation of the Einstein tensor gives 

2<5Goo = 2zh k k - 4:e- 2z k 2 zh 01 + e - 2z k 2 (h kk - fin), 

28Gq\ = ikhn — ikh kk + ik(4z + 6z 2 )/iqi, 
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2SGoa = ikhia + ik(4:Z + 6z 2 )h aa — ike 2z k 2 ho a , 
25Gn = e 2z [h n + 3zh u - (42 + 6z 2 )h n - h kk - 3zh kk ], 
25G\ a = e 2z [hi a + 3zh\ a — (4z + 6z 2 )hi a ] — k 2 h$ a — k 2 zho a , 
2SG 22 = e 2z [h 22 + 3zh 22 - (4z + 6z 2 )h 22 - h kk - 3zh kk ] + 2k 2 h 01 + 2k 2 zh 01 + k 2 (h n + h 22 - h kk ), 
26G 33 = e 2z [h 33 + 3zh 33 - (4z + 6z 2 )^ 33 - h kk - 3zh kk ] + 2k 2 h m + 2k 2 zh m + k 2 (h n + h 33 -h kk ), 

2SG 23 = e 2z [h 23 + 3zh 23 - (4z + 6z 2 )h 23 ] + k 2 h 23 ; 
and variation of the equation of motion of the continuum gives 

(e + <t)Aoi = {2n + 3A + 5a)zh 01 - (fj, + a)hu - i(A + a)h kk , 

(e + a)h 0a = (2[i + 3A + 5<r)zh 0a — (/i + a)h la . 

The index a, appearing in the expressions for 8G^ U as well as in the varied equations of motion 
of the continuum, assumes values 2 and 3. The expressions for 5T^ U and 25G fJiL , are given modulo 
factor e %kxl ; this is, however, of no significance since we aim to equate them anyway. 

When writing down the varied Einstein equations we can use formulas for the background 
solution to eliminate both po and z from our equations (and keep z, z, e+a, A + er and /i + er only). 
Again, as in the background solution, we find that some terms of the varied Einstein equations 
are just time derivatives of other equations. In addition to the varied equations of motion of the 
continuum cited above, we obtain the following set of independent equations: 

2zh kk + e~ 2z k 2 (h kk - h xl - 4zh i) = -^e~ 3z (e + <r)h kk , 

hn - h kk = e~ 3z (e + cr)h 01 , 
h la - e~ 2z k 2 h 0a = e~ 3z (e + <r)h 0a , 
(h 22 - h 33 ) + 3z(h 22 - h 33 ) + e- 2z k 2 (h 22 - h 33 ) = -e- 3z (^ + n){h 22 - h 33 ), 
h 23 + 3zh 23 + e- 2z k 2 h 23 = -e- 3z {n + cr)h 23 . 

3.3 Modes 

The set of equations we have found is evidently not totally coupled; it can be rather divided 
into several subsets of coupled equations. Physically, this means that there are different modes 
of wave propagation for the given wave covector. Since general relativity without matter and 
classical elasticity are both limit cases of relasticity, we expect that there will be two modes of 
gravitational waves, two modes of transversal acoustic waves and one mode of longitudal acoustic 
waves. Also, there should remain one odd mode coming from the freedom of choosing the zero 
time hypersurface. 

First we have equations for h 22 — h 33 and h 23 , which can be written in a compact form as 

hjj + 3zhl + [e- 2z k 2 + e- 3z (fi + a)] h% = 0, (5) 
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where tiL is the tensor part of (a symmetric 3x3 matrix with nonzero components h 22 = 
—h 33 = (/i22 — h 33 )/2 and h 23 = h% 3 ). These equations clearly describe gravitational waves with 
polarizations © (the component h 22 ) and ® (the component /1I3); or, speaking in cosmological 
terms, the tensor perturbations (see |6j)- The equations coincide with those valid in general 
relativity for an universe filled with an ideal fluid, see [TO] , if /i + a vanishes. To see where this 
condition comes from, note that for an ideal fluid the energy per particle e depends only on 
the particle density n. If we write n = V detH AB , we can compute the second derivatives of 
e = e(n) with respect to H AB and compare them with the expression in subsection 13.11 In this 
way we obtain /i = —nde/dn = (l/3)de/dz. However, the identity for de/dz, to be found also in 
subsection 13.14 implies that the right hand side equals —a. Thus, an ideal fluid can be defined 
as an elastic continuum with /1 + a = 0. 

Next we have two pairs of equations containing hg a and h\ a only, 

• 2^ + 3A + 5cr fi -\- a 

"Oa — ; ZtlQa ; u,\ai \0) 

e + a e + a 

h la = [e~ 2z k 2 + e~ 3z (e + a)] h Qa . (7) 

Consider equations for h,Q2 and hyi. Since the perturbation of g\i represents the shear strain 
in the x 2 direction and the perturbation of (702 represents the motion of the particles of the 
continuum in the same direction (for an observer for whom the time t is locally synchronized), 
this mode can be identified with the transversal acoustic wave that oscillates in the x 2 direction. 
Analogically, equations for /i 3 and /113 describe the propagation of a transversal acoustic wave 
that oscillates in the x 3 direction. These modes belong to the sector of vector cosmological 
perturbations. 

There are three independent combinations of components of h left, namely /ioi, ^11 and hkk, 
describing scalar cosmological perturbations. But they are all inflicted by the remaining freedom 
of choosing the zero time hypersurface. Pure shift of this hypersurface, corresponding to a plane 
wave with the wave covector k — {k, 0,0), yields a contribution to h proportional to hoi = 1, 
h\x = h 2 2 = h 33 — 2z. To identify physically relevant perturbations, we define combinations 
of components of h that are invariant with respect to this hypersurface change. We do this by 
introducing yoij Uii an d y such that 

hoi = yoi + y, 
hn = yu + 2zy, 
hkk = yu + Gzy. 

The interpretation of these formulas is as follows. If z ^ 0, we can always choose the zero time 
hypersurface in such a way that hkk — hn. For such hypersurface, yoi and yu are equal to hoi 
and hu respectively. If we want to choose another hypersurface, we simply add the hypersurface 
shifting terms proportional to y. Since we still want t to be the proper time of particles, we will 
have differential equation for y, but we can see that the physical information about the longitudal 
acoustic waves is contained in yoi and yu only. 

Using the definition of y's, we obtain a set of equations that can be neaten to one equation 
for y, 

y=-^-^{e + cj)yoi, (8) 



8 



and a pair of equations not containing y which invariantly describe the propagation of a longitudal 
acoustic wave, 



Vox 



2^ + 3A + 5er. 1 _ 3z 



z + —e~ 6z {e + o) 



Vn 



e + a 4i 



1 A + 2fi + 3a 

yoi — o — 7T yn > ( 9 ) 

z e + a 



yoi-^e 32 (£ + a)t/ n . (10) 



4 Example 

Consider a continuum with the stiff ultrarigid state function introduced in [9]. The authors 
define the principal linear particle densities ni, n 2 , n 3 in terms of H AB , and write the state 
function as 

p = A [(nm 2 ) 2 + (n 2 n 3 ) 2 + (n 3 m) 2 ] + B, 

where A and B are constants. If we choose, for simplicity, A ~ 1 and B — 0, the state function 
as a function of H AB is 

g =l(H cc H DD -H CD H CD ). (11) 

If we treat all H AB, s as truly independent variables, the second term depends on the antisym- 
metric part of H AB ; however, there is no need to modify it as long as we use our formalism of 
symmetrized higher derivatives. 

The first and second derivatives of g with respect to H AB are 

®6 r TJ CC x x zjCD 



Q H AB 

d 2 g 



= 6 AB H°° - S ac S B dH 01J , 
— SabScd — 8ac$bd- 



dH AB dH CD 

Evaluating g and its derivatives for H AB = e~ 2z 6 AB we get 

_ 4z dg _ 2 d 2 g 

6= ' dH AB= AB ' dH AB dH CD = bABdcD ~ b AC(>BD, 

and, consequently, 

p = 0, e = 3e~ z , a = e~ z , \ = -e~ z , p = e~ z . (12) 

The equation 6i 2 = po + e~ 3z e reads 

and has an (expanding) solution 

z = -\n2 + -\nt. (13) 
4 2 

This is just the well known formula a oc \ft for a universe filled with an ideal fluid with the 
equation of state a = (l/3)e (or, in more common notations, p — (l/3)p). 

Using z = (2i) -1 and e~ 2z — (V^t)^ 1 , we can express coefficients in the equation, or system 
of equations, for each kind of waves as functions of time. The equation for gravitational waves 
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reads 



(14) 



Besides of that, we have a pair of equations for transversal acoustic waves 




( 



k 2 
V2t 



+ 




and a triplet of equations for longitudal acoustic waves, 



2 

v = --yoi, 



3 

2/01 = 7772/01 





The first-order equations for h$ a and /ii Q , as well as those for yoi and yu, can be combined into 
a second-order equation for either variable. We choose the variables ho a and yoi, because the 
corresponding equations contain k 2 only as a factor in front of the variable itself, thus becoming 
hyperbolic partial differential equations when one passes from k 2 back to the Laplacian. The 
equations read 



A remarkable feature of these equations is that they coincide in the limit of small k. Thus, if we 
use for the description of longitudinal and transversal acoustic waves the functions ho a and yoi , 
the waves behave identically in the regime in which they stay well outside the horizon. There is no 
point in extending the comparison to gravitational waves, since the field tif^ in the corresponding 
equation is of completely different nature than the fields ho a and yoi in the equations for acoustic 
waves (which are at least akin). Note, however, that if we pick the variables hi a and yu instead 
of hoa and yoi to describe acoustic waves, or multiply ho a and yoi by e~ 2z (the same factor that 
appears in the definition of h%j in terms of Sgij), the resulting equations will coincide not only 
with each other, but also with the equation for hfj in the limit of small k. 

For both gravitational and acoustic waves, oscillating fields as functions of time are depicted 
in figure [U The curves denoted by g ('gravitational'), t ('transversal') and I ('longitudinal') are 
graphs of the functions h 22 or h 23 , hoi or ho3, and yoi respectively. For comparison, graphs of the 
functions h 22 or h 23 and yoi for an ideal fluid with the equation of state a = (l/3)e are included 
into the figure, too (denoted by g fluid and I fluid', of course, there are no transversal acoustic 
waves in an ideal fluid). All curves are computed for k = 1, with the oscillating field equal to 
1 and its derivative vanishing at the moment t — 0.01. The waves cross the horizon at t ~ 1; 
before that their wavelength is larger and afterwards it is smaller than the horizon length. The 
main new effect of nonzero shear stresses, distinguishing the continuum under consideration from 
an ideal fluid, are oscillations of the fields at superhorizon scales. Asymptotically, the oscillating 
factor is cos[(v / 15/4) logt + tp] for all three kinds of waves. 




(15) 




(16) 
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Figure 1 : Time evolution of gravitational and acoustic waves 

5 Conclusion 

We have developed a formalism describing the propagation of plane waves in a general rel- 
ativistic homogeneous and isotropic flat continuum. In comparison with [5], where the same 
problem is studied, we did not include into the theory small perturbations to the Euclidean 
metric of the material manifold (tensor b a b in the notation of [5]). In this way we simplified the 
matters as much as possible. 

We were successful in identifying relativistic generalizations of acoustic waves. The longitu- 
dinal mode was encumbered with the 'time wave' gauge freedom, and the definition of gauge 
invariant y's was needed. After introducing them we were able to describe the wave by two 
coupled first-order differential equations. This is equivalent to one second-order equation, hence 
our theory is a straightforward generalization of its nonrelativistic counterpart. Again, this can 
be compared with [5] where the propagation of the longitudinal wave is described by two coupled 
second- order equations (equations (3.12) and (3.15) in subsection 3.1). The cited work differs 
from ours in that respect that it uses Newtonian gauge. Of course, the transition from one gauge 
to another cannot increase the number of degrees of freedom, therefore in Newtonian gauge 
the longitudinal wave must be also described by two first-order equations (imposing constraints 
on initial conditions for the second-order equations). Perhaps the comoving proper-time gauge 
we have chosen for our calculations is better adapted for derivation of these more fundamental 
equations. 

We kept the formalism as easy and simple as possible, although sometimes we have preferred 
comprehensibility to simplicity. For example, the introduction of po and e was not necessary. We 
could use g (and scalars entering the derivatives of g with respect to H AB ) only, but we felt that 
the quantities e (and a, /i and A) were more appropriate for the description of the continuum. 
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Also, we introduced z and h's just to simplify the expressions even if the meaning of hoi is far 
from obvious. And our choice of y's reflected the most natural initial condition h\\ = hkk- On 
the other hand, the fact that the resulting equations do have singularity for z = suggests that 
this condition is not the best. For those interested in nonsingular equations, note that using 

hox = v i + v, 

hkk = v k k + 6zv, 
Ae~ 2z k 2 + 2e- 3z (e + a) 

^11 = ^2 V ^ + 2ZV > 



we obtain equations 



. _ e- 3z (e + 2^ + 3A + 6o-) 

6 e - 3z ( £ + 2 M + 3A + 6a) . _ 2 2 

Vkk = , _ 2z]2 , q , z( — ; — rzvkk + 2e k v m , 



2[i + 3A + 5er 
u i = ; zvoi 



12/z + A + 3ct 2e~ 6z (n + a) e - iz (e + + 3A + 6a 



Vkk- 



e + a Ae- 2z k 2 4e- 2z k 2 + 3e- 3z (e + a) _ 

However, these equations look no simpler than the equations for h's, and we used z ^ when 
deriving them anyway. 

When applying our theory to a stiff ultrarigid continuum, we found that the time evolution of 
longitudinal and transversal acoustic waves, if described by the functions ho a and yoi, coincide 
in the limit of small k. This is not a unique feature of the state function considered. The same 
statement holds for any 'linear' continuum, that is, continuum with po = and a, A and /i 
proportional to e, provided the constraint on the coefficients of proportionality following from 
the relation da/dz = —(2a + 3A + 2/x) is satisfied. 
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